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ABSTRACT 



If the angular momentum of the molecular cloud core were conserved during 
the star formation process, a new-born star would rotate much faster than its 
fission speed. This constitutes the angular momentum problem of new-born 
stars. In this paper, the angular momentum transfer in the contraction of a 
rotating magnetized cloud is studied with axisymmetric MHD simulations. 
Owing to the large dynamic range covered by the nested-grid method, the 
structure of the cloud in the range from 10 AU to 0.1 pc is explored. First, 
the cloud experiences a run-away collapse, and a disk forms perpendicularly 
to the magnetic field, in which the central density increases greatly in a finite 
time-scale. In this phase, the specific angular momentum j of the disk decreases 
to ~ 1/3 of the initial cloud. After the central density of the disk exceeds 
~ 10 10 cm~ 3 , the infall on to the central object develops. In this accretion stage, 
the rotation motion and thus the toroidal magnetic field drive the outflow. The 
angular momentum of the central object is transferred efficiently by the outflow 
as well as the effect of the magnetic stress. In 7000 yr from the core formation, 
the specific angular momentum of the central 0.17M© decreases a factor of 10~ 4 
from the initial value (i.e. from 10 20 cm 2 s _1 to 10 16 cm 2 s _1 ). 



Subject headings: ISM: clouds - ISM: jets and outflows - ISM: magnetic fields - 
stars: formation - stars: rotation 
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1. INTRODUCTION 

The problem of star formation has been the object of serious study for many 
years. In the star formation process, almost all the angular momentum should be 
removed from the molecular cloud in the course of forming stars. For example, using 
rotation periods and the radius of classical T Tauri stars (CTTS) [P = 3 — 10 days 
and i?* ~ 2R Q (Bouvier et al. 1993)], the specific angular momentum of such stars 
is approximately equal to j ~ AnR 2 J(5P) ~ 5.6 x 10 16 cm 2 s _1 ( J R,/2 J R o ) 2 (P/10day)- 1 , 
where we assumed the star rotates uniformly. On the other hand, if the velocity 
gradient of 0.3km s _1 pc _1 — 4km s _1 pc _1 (Goodman et al. 1993) observed in 
NH 3 cores (density ^ 10 4 cm~ 3 ; size ~ 0.1 pc) comes from the rotation, j reaches 
~ 5 x 10 21 cm 2 s 1 (i?/0.1pc) 2 (f2/4km s _1 pc -1 ). Therefore, the specific angular momentum 
should decrease a factor of 10~ 5 when a star is formed from a molecular cloud core. 
Assuming j is conserved through the collapse, the gravity is balanced with the centrifugal 
force at the centrifugal radius as R c = j 2 /{GM) ~ 0.06pc(j/5 x 10 21 cm 2 s _1 ) 2 (M/lM )- 1 . 
This leads to the conclusion that if the angular momentum is not efficiently reduced or 
transferred, the cloud could not shrink to become a stellar object (R c 3> R*). This is called 
the angular momentum problem for new-born stars. 

To reduce the spin angular momentum of newly formed stars, two mechanisms are 
proposed. One is to transfer the cloud spin angular momentum to the orbital angular 
momentum of binary or multiple protostars if they are made by fission. The other possibility 
is magnetic braking, in which angular momentum is transferred to the ambient medium 
by the magnetic stress (Spitzer 1978). The latter works even for the case of a single star. 
In the case of the parallel rotator, in which the directions of magnetic fields and angular 
momentum coincide with each other, the time scale of deceleration of the angular rotation 
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speed (Q) is given by 



^ ° 2"7rGa 1/2 



where a, p a , Va, and G denote, respectively, the column density of the molecular core, the 
ambient density, the Alfven speed of the ambient medium, and the gravitational constant 
(Ebert, von Hoerner, & Temesvary 1960, Mouschovias, & Paleologou 1980). As long as the 
molecular core evolves in a quasistatic manner, the decrease in the angular momentum L 
is well fitted by the above time scale as L = L(t = 0) exp(—t/t B ) (Tomisaka, Ikeuchi, & 
Nakamura 1990). Since the first factor 2nGa/Bo is approximately equal to or larger than 
unity for a gravitationally contracting cloud (Tomisaka 1995), ts is longer than the free-fall 
time of the ambient low-density medium. This raises the question of whether or not the 
angular momentum is efficiently transferred while the contraction proceeds. 

Recently, bipolar outflows have been found in a wide variety of objects which indicate 
star formation (for a review, see Bachiller 1996). Based on the magnetic acceleration model 
of outflow (Blandford, & Payne 1982; see also Konigl 1989 and Wardle, & Konigl 1993), 
magnetic fields extract the angular momentum from the disk and transfer it to the outflow 
gas. Therefore, the outflow gas is able to carry the excess angular momentum directly. 

From MHD simulations of the gravitational contraction of a magnetized, rotating 
isothermal molecular cloud, it is shown that the outflow begins after the central density 
exceeds 10 10 cm -3 and an adiabatic core forms in the center (Tomisaka 1998, hereafter 
Paper I). The outflow is never driven in the run-away collapse phase, since rotation motion 
and thus toroidal magnetic fields are unimportant in this phase. Therefore, in this paper, 
we study the evolution of the cloud from the runaway collapse phase to the accretion 
phase, continuously. In particular, we try to determine when the redistribution of angular 
momentum occurs and whether the outflow plays an important role in extracting the 
angular momentum from the contracting gas. The model and the numerical method are 
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described in §2. In § 3, the evolution of the specific angular momentum distribution is 
shown. In addition, a comparison is made between the effects of the magnetic torque and 
outflow to reduce the angular momentum. We discuss the ideal MHD condition and the 
evolution of late stage in §4. 



Model and Numerical Method 



We began our simulation from an infinitely long, cylindrical, rotating, isothermal 
cloud in hydrostatic balance. In terms of the gravitational potential, ipo, the magnetic flux 
density, B = B z e z , the angular rotation speed, Q , and the density p , the hydrostatic 
configuration is governed by the equations 



rQ 2 



# c 2 dp 1 dB 2 = Q 



dr 
1 d 



Po dr 
di/j \ 



8ir dr 
= AtcGpq. 



(2) 



r dr \ dr 

These equations have a solution for the radial distributions of density, rotational velocity, 
and magnetic flux density as follows (Stodolkiewicz 1963): 

-2 



(3) 



Po(r) = Pc 



r 2 /(8H' 



-1/2 



(4) 

(5) 
(6) 



V<t> (r) = rQ {r) = rQ c [l + r 2 /{8H 2 ) 

B z (r) = B c [l+r7(8# 2 )p, 

where H represents the scale-height as AnGp c H 2 = c 2 + B 2 /(87ip c ) + 2VL 2 H 2 [for 
Vt c < (2irGp c ) 1 ^ 2 ]. Here, quantities with a subscript c denote those for the cloud 
center (r = 0). This hydrostatic solution contains two nondimensional parameters: 
to = n c /(47vGp c y/ 2 and a = B 2 /(4iTp c c 2 ). 

We assume that the gas obeys the ideal MHD equations. The gas is assumed isothermal 
for low densities. However, after the density exceeds ~ 10 10 cm~ 3 , the interstellar gas is 
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not efficiently cooled and behaves adiabatically (Larson 1969). To mimic this situation, a 
barotropic relation such as 



PS 



= c 2 s p, . . . (p < p crit ) ^ 

= C 2 s p crit (p/pcrit) r , • • • (P > Pcrit) 



was adopted without using the energy equation, where we take p crit = 10 10 H 2 cm -3 . The 
polytropic exponent T was taken as 5/3 and 2. (The specific heat ratio 7 never exceeds 5/3. 
However, we found that the flow outside the core and the mass-accretion rate are essentially 
the same irrespective of the values of the core's polytropic exponent, T. Therefore, we 
used the larger exponent to reduce the effort to perform numerical computations.) To 
initiate gravitational contraction, we added density perturbation with small amplitude. 
The periodic boundary condition was applied to the boundaries in the z-direction, and 
the separation between two boundaries was chosen to be identical with the wavelength 
of the gravitationally most unstable mode (Matsumoto, Hanawa, & Nakamura 1997). To 
ensure the fine spatial resolution necessary especially near the forming adiabatic core, the 
nested-grid method was applied to the MHD finite difference scheme (Tomisaka 1996a, 
1996b, 1998). In the nested-grid method, 15 levels of grids are used from LO (the coarsest) 
to L14 (the finest). A grid spacing of Ln was chosen as 1/2 of that of Ln — 1, and each 
level of grids was divided into 64 x 64 cells. The code was tested by comparing the result of 
gravitational collapse calculated by a code without the nested-grid method (for details, see 
Tomisaka 1996b). 



3. RESULTS 

Here, the results of the model with w — 1/2 and a = 1 are shown. This corresponds 
to Q — 2.78n^4 2 km s _1 pc _1 and B c = 13.3n^4 2 c s i9o pG. Here, n c4 and c s i 90 represent 
respectively p c /10 4 H 2 cm~ 3 and c s /190m s" 1 . The parameters are taken to be the same as 
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the model shown in Paper I. Evolutions are also the same. 



3.1. Specific Angular Momentum Distribution 



In Figure 1, the specific angular momentum j is plotted against the mass measured 
from the center. Since the density increases monotonically reaching the center, a mass shell 
with higher density is considered to be located nearer to the center than that with lower 
density. The total mass of the gas whose density is higher than p 1: 



are calculated at three epochs. Figure 1 is made by plotting j(> pi) = L/M against 
M(> p\) for various p\. This is similar to the specific angular momentum spectrum (e.g. 
Fig. 5 of Norman, Wilson, & Barton 1980) in which M(< j), the total mass whose specific 
angular momentum is smaller than j is plotted against j. However, the plot in Figure 1 has 
an advantage that the gas with smaller M is necessarily located nearer to the center. 

Open circles represent the j{< M) distribution in the early phase [ 
t = 5.214(4vrG'p c )- 1 / 2 ~ 0.9n; 4 1/2 Myr ], in which the central density increases a 
factor 10 from the initial state (i.e., p c reaches ~ 10 5 n c4 H 2 cm -3 ). At this stage, the density 
distribution is almost spherical. Filled squares denote the distribution at the end of the 
run-away collapse phase (t < 5.980(47rG'p c ) _1 / 2 ~ 1.05n~4^ 2 Myr; the structure at the stage 
is shown in Fig. 2a of Paper I). Since the motion crossing the magnetic fields is blocked, a 
disk running perpendicularly to the magnetic fields is formed. From Figure 1, it is shown 
that the angular momentum contained in the central ~ 0.5c 3 190 n~ 4 1//2 M Q is reduced to 
~ 1/3 before a core is formed. 




(8) 



and the correspondent total angular momentum, 




(9) 
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This reduction is not owing to the angular momentum transfer by the magnetic fields 
but to the change of the shape of high-density region (from spherical to disk-like). This 
is confirmed by the fact that the specific angular momentum spectrum (not shown in this 
Letter) is not changed during the run-away collapse phase. The disk is mainly formed by 
the gas flow parallel to the magnetic fields (z-direction). Therefore, central part of the 
run-away collapsing disk is made by a gas with small angular momentum sitting initially 
near the z-axis. In other words, to the high-density region the mass is preferentially 
gathered rather than the angular momentum. The difference between open circles and 
filled squares in Figure 1 is owing to the segregation of low-angular momentum gas in the 
run-away collapse stage. However, the amount of the angular momentum segregation is too 
insufficient to explain the angular momentum problem. 

The specific angular momentum, j, is approximately proportional to the accumulated 
mass, M. This relation coincides with the prediction for the thin disk by Basu(1998). This 
relation is explained from the solution of the contracting isothermal thin disk, i.e., the 
column density a ~ c s /Gr and the angular rotation speed Q oc r~ x (Matsumoto et al. 1997) 
lead to the specific angular momentum j = Qr 2 oc r oc M = / ardr. 

After the central density of the disk exceeds p cr it ~ 10 10 H2cm~ 3 , a new phase appears: 
the accretion phase. Thermal photons emitted from the dust which cooled the cloud become 
optically thick; as a result, an almost spherical adiabatic core is formed (in the adiabatic 
core, the thermal pressure becomes more important). Gas around the core begins to accrete 
on the core. 

As shown in Paper I, after ~ 10 3 yr has passed,^ the outflow begins (Figures 
2b and 3 of Paper I). The size of the seeding region (origin of the outflow) is 

2 This age is independent on the values of the assumed polytropic index, T = 5/3 and 
2. Furthermore, the mass outflow rate differs slightly: T = 2 model shows approximately 20 
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r ~ 5 x 1CT 3 c s I (AnGpc) 1 / 2 ~ 35c s i 9 on~ 4 1//2 AlJ. The seeding region expands radially outward. 
This outflow is driven by the gradient of the magnetic pressure of the toroidal magnetic 
fields, — V-B|/87r, which are made by the rotation motion (Paper I;Kudoh, Matsumoto, & 
Shibata 1998). The magnetic fields exert torque on the outflowing gas to increase its angular 
momentum. On the other hand, they exert torque on the disk to decrease the angular 
momentum (Blandford, & Payne 1982). In r = 0.04(4ttG'p c )~ 1/2 ~ 7000< 4 1/2 yr (filled 
triangles; r is the time elapsed from the core formation epoch), the outflow expands and 
reaches z ~ 0.25c s (47rGp c ) _1 / 2 ~ 1800c sl90 n C 4 /2 AU. The angular momentum distribution 
at that time is shown by the solid line in Figure I. It is shown that the specific angular 
momentum contained in the central ~ O.17c^ 19O n C 4^ 2 M (p > 10 10 n C 4H 2 cm~ 3 ) has been 
reduced to a factor ^ 1CT 4 of the initial value. 

3.2. Angular Momentum Flux and Magnetic Torques 

Consider the total angular momentum, L(< M), contained in the mass M. The 
increase and decrease of the angular momentum are caused by (1) the inflow (accretion 
flow) mainly in the disk dL/dt[ in) (2) the torque exerted on the gas disk from the magnetic 
fields N, and (3) the outflow from the disk dL/dt\ out : dL/dt — dL/dt\- m + dL/dt\ out + N. It 
is to be noted that only the first term is positive. To see which term is important to reduce 

% smaller outflow rate than the model of T = 5/3. However, in the case of slowly rotating 
cloud as a = 1 and uo = 1/10, which leads to a smaller centrifugal radius, in a model 
assuming T = 2 the outflow begins earlier than than that of T = 5/3. This seems to indicate 
that unless the centrifugal radius from which the outflow is mainly ejected is formed near 
the center, especially, in the hard core, assuming such a relatively hard polytropic index is 
justified. Quantitative comparison will be made in a forthcoming paper. 
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the angular momentum, we plot these three quantities for the three epochs in Figure 2. To 
calculate dL/dt[ m and dL/dt\ out , we consider a cylinder which covers the isodensity surface. 
The dL/dt\ OVLt is measured by the angular momentum flux convected outwardly through the 
upper and lower surfaces of the cylinder, and dL/dt[ m is the flux running radially inward 
through the side surface. 

Before the core formation (circles and squares), the angular momentum inflow 
(dL/dt\ in ) is larger than the magnetic torque (N) which reduces the angular momentum. 
(This does not mean that the specific angular momentum increases. Note that dM/dt m is 
also positive.) In contrast, in the accretion phase after outflow blows, the magnetic torque 
(filled triangles) is comparable with dL/dt\- m (open triangles) near M ^ 0.17c 3 sl90 n^ 2 M Q 

I/O 

(r = 7000n c4 yr). Further, from the figure, it is shown that the angular momentum 
transport by outflow dL/dt\ out (open stars) is more efficient than that of the magnetic 
torque. By these two effects, the angular momentum transferred by the inflow inside the 
disk is totally removed for central ~ 0.17c 3 igo n~^ 2 M Q . As a result, j(< 0. 17 'c 3 igo n~l^ 2 M & ) 
is reduced a factor of 10 -4 from the initial value. 



4. DISCUSSION 

The size of the outflow reaches z ~ 0.25c s /(47rGp c ) 1 / 2 ~ 1800c s i 90 n~y 2 A\] in 
t ~ 7000n~4^ 2 yr. However, CO observations show that the molecular outflows have typical 
spatial size of 8000 AU - 4 pc (Fukui et al. 1993). This simulation is limited in time and 
should be extened further. However, calculation becomes harder after the core mass exceeds 
~ 0.1c 3 l90 n c l^ 2 M Q due to high density at the cloud center. Let us consider whether the 

— 1/2 

magnetic angular momentum transport is effective still in the later stage r > 7000n c4 yr. 
As shown in Figure 1, j is approximately proportional to M at the epoch of the 
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core formation. Writing down the relation as j = p x GM/c s , the coefficient p is 
approximately equal to O.f for the model shown in Figure 1. If j is conserved, the 
centrifugal radius (R c ), at which the gravity is balanced with centrifugal force, is given by 
R c = f/GM = p 2 GM/c 2 s ~ 25O(p/O.l) 2 (M/lM )c; 1 2 9O AU [a similar relation is pointed out 
by Basu (1998)]. This indicates that the centrifugal radius expands with time. On the other 
hand, a self-similar solution for the rotating isothermal thin disk has been found by Saigo & 
Hanawa (1998). Their solution in the accretion phase also shows that the radius where 
takes the maximum increases with time as oc c s t. Since the outflow is ejected from a radius 
where is important (Paper I), these two indicate that the seeding region of the outflow 
moves outwardly. Using the solution by Saigo & Hanawa (1998; for an example, see their 
Figure 7), the surface density (£) and thus the density (p) decrease as collapse proceeds 
and the seeding region moves outwards. This shows that the density near the seeding region 
decreases as the collapse proceeds. 

This leads to the conclusion that the coupling between gas and magnetic fields (Nakano 
1990) becomes stronger as long as we consider the seeding region, indicating that the 
mechanism of angular momentum transfer works also in the later stage of the evolution. 

The author thanks R. Ouyed, the referee, for his fruitful comments to improve 
the paper. He also thanks T. Hanawa, T. Matsumoto, and F. Nakamura for useful 
discussion. This work was partially supported by Grants-in-Aid from the Ministry of 
Education, Science, Culture, and Sports (10147105, 11640231). Numerical calculations were 
performed by Fujitsu VPP300/16 at the Astronomical Data Analysis Center, the National 
Astronomical Observatory, and NEC SX 4/2B at the Integrated Information Processing 
Center, Niigata University. 
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Fig. 1. — Specific angular momentum j is plotted against the accumulated mass from the 
center M. Open circles, filled squares, and triangles denote, respectively, the beginning and 
the end of the run-away collapse phase and the accretion phase. 
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Fig. 2. — The torque exerted on the cloud core by the magnetic fields, N, the angular 
momentum inflow rate convected by the inflow in the disk, dL/dt\ in , and its outflow rate 
transferred by the outflow, dL/dr\ ont . Filled and open symbols denote respectively the 
distribution of N and that of dL/dt\- m . The circles, and squares denote, respectively, the 
beginning and the end of the run-away collapse phase. Open stars represent the angular 
momentum outflow rate dL/dt\ out in the accretion phase. 



